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Abstract. This paper directly builds upon previous work in [7], where we introduced new 
reduced basis a posteriori error bounds for parametrized saddle point problems based on Brezzi's 
theory. We here sharpen these estimates for the special case of a symmetric problem. Numerical 
results provide a direct comparison with former approaches and quantify the superiority of the 
new developed error bounds in practice: Effectivities now decrease significantly; consequently, the 
proposed methods provide accurate reduced basis approximations at much less computational cost. 
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Introduction. The reduced basis metliod is a model order reduction approach 
that permits the cfhcient yet reliable approximation of input-output relationships in- 
duced by parametrized partial differential equations. In contrast to generic discretiza- 
tion techniques where approximation spaces are not correlated to the physical proper- 
ties of the system, the method recognizes that the solutions to a parametrized partial 
differential equation (PDE) are not arbitrary members of the infinite-dimensional 
solution space but rather reside or evolve on a much lower-dimensional manifold. 
Exploitation of this low-dimensionality is the key idea of the reduced basis approach. 

Designed for the real-time and many-query context of parameter estimation, op- 
timization, and control, the method provides rapidly convergent and computationally 
efficient approximations equipped with practicable and rigorous error bounds. Built 
upon a high-fidelity "truth" finite element discretization, the reduced basis approxima- 
tion is defined as a Galerkin projection onto a low-dimensional subspace that focuses 
on the solution manifold induced by the parametrized PDE. The error in the reduced 
basis approximation is then measured relative to the "truth" problem formulation 
and can be quantified by rigorous and inexpensive a posteriori error bounds. To the 
method's key features also belongs an Offline-Online computational strategy that en- 
ables the highly efficient (Online) computation of both reduced basis approximations 
and error bounds for any parameter query at the expense of increased pre-processing 
(Offline) cost. Finally, reduced basis approximations and error bounds are intimately 
linked through a greedy sampling approach, in which the (Online-)inexpensive error 
bounds are used to construct the reduced basis approximation spaces more optimally. 

This paper directly builds upon our work in 7 , where we introduced new re- 
duced basis a posteriori error bounds for parametrized saddle point problems based 
on Brezzi's theory [3l|4]. In contrast to former approaches based on Babuska's the- 
ory for noncoercive problems [U 1121 , these do not involve the highly expensive 
estimation of the Babuska inf-sup stability constants but only much less expensive 
calculations; as separate upper bounds A" and A^* for the errors in the reduced basis 
approximations for the primal variable u and the Lagrange multiplier p, respectively, 
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they moreover enable the systematic estimation of engineering outputs depending on 
either of the two. Numerical results showed that the bounds A" were reasonably 
sharp and thus very meaningful. However, overestimated the actual error in the 
reduced basis approximations for p rather pessimistically; here, error estimates A^'^ 
based on Babuska's theory achieved much better results. 

In this paper, we now investigate symmetric parametrized saddle point problems. 
We find that in this special case, a posteriori error bounds A", A^ developed in [7] 
may be improved considerably. Considering the same Stokes model problem as before, 
numerical results quantify this behavior in practice through a direct comparison with 
former techniques: Effectivities for both A" and A'' now decrease significantly, and 
perform in neither case worse than effectivities associated with A^'^. 

The paper is organized as follows: In [JTJ we recall the general formulation of 
a symmetric parametrized saddle point problem together with its "truth" approx- 
imation. Section [5] then describes the reduced basis method with a strong focus 
on specifics to the symmetric case considered. In ^2.\\ we define the reduced basis 
approximation as the Galerkin projection onto a low-dimensional reduced basis ap- 
proximation space. A priori as well as a posteriori error estimates shall be specialized 
to the symmetric case in i i2.2l Reduced basis approximations and a posteriori error 
bounds can be computed (Online-)e5iciently as summarized in i i2.3l This enables 
us to engage adaptive sampling processes for constructing computationally efficient 
reduced basis approximation spaces, which shall be outlined in i j2.4l In [J3l numerical 
results demonstrate the performance of the improved a posteriori error bounds in 
practice. Finally, in [21 we give some concluding remarks. 

1. General Problem Statement. 

1.1. Formulation. We shall briefly recall the setting introduced in [7]. Let 
Xf. and Yc be two Hilbert spaces with inner products (•, ■)xe, {'j')y^ and associated 
norms \\-\\x^ = \l (•, ■)x^, II • lli'c = V (•, ■)y^-, respectively^ We define the product space 
= XeXFe, with inner product {■,-)z^ = (•, Oxe + (•, O^e and norm ||-|j2^ = ^J{-.^ ■)z,. 
The associated dual spaces are denoted by X'^, FJ, and Z'^. 

Furthermore, let V C M" be a prescribed n-dimensional, compact parameter set. 
For any parameter /i G I?, we then consider the continuous bilinear forms a(-, •;/!) : 
XcX Xc^R and b{-, •; n) -.XcXY^^ mH 

7g(A*) = sup sup ||"||^'n'^^ < oo, VfieT), (1.1) 

uex^vex^ \\u\\xjv\\x, 

^lin) ^ sup sup /i^'i;'^^ < oo, V /i e P. (1.2) 
We moreover assume that a(-, •; /i) is coercive on Xc, 

II V II 

and that b{-, •; /i) satisfies the inf-sup condition 

EE inf sup /.^'^'ll'l'^ > 0, V A* e P. (1.4) 



<(/i) = inf ,',,2 > 0, V ^ e P, (1.3) 



^Here and in the following, the subscript e denotes "exact". 

^For clarity of exposition, we suppress the obvious requirement of nonzero elements in the de- 
nominators. 
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We now consider the following variational problem: For any given yu G 2?, we find 

(uo(/i),Pc(M)) G -''^c ^ Yc such that 

b{uc{fi),q;fj.) ^ g{q;fi), VgeFc, 

where f{-',fJ') and g{-', IJ-) are bounded linear functionals in and Y^, respectively. 
From the results of Brezzi |3l[4], it is well-known that under the assumptions (jl.ip . 
(|1.2p . (|1.3p . and (|1.4p . the above problem (jl.Sp is well-posed and has a unique solution 
for any f e X^, g{-,fi)eY^. 

In contrast to this very general setting considered in [7], we here additionally 
assume that the bilinear form a(-, ■; ^) is symmetric for any fj. € V. As a continuous, 
symmetric, and coercive bilinear form, a(-, •;/!) now defines an inner product on Xg 
for any ^ € V with an associated norm \\ ■ \\x^,fj, = -\/a(-, •; /x) equivalent to || • \\x^- 
We note (see, e.g., [HH]) that the solution {uc{^J,),Pe{n)) to (|1.5p then corresponds to 
a saddle point of the Lagrangian functional 

C{v, q; fi) = la{v, v; fi) + b{v, q; n) - f{v; fi) - g{q; fi), V {v, q) £ Z^. 

1.2. Truth Approximation. We now introduce a high-fidelity "truth" approx- 
imation upon which our reduced basis approximation will subsequently be built. To 
this end, let X and Y denote finite-dimensional subspaces of Xc and Yc, respectively. 
We define the product space Z = X x Y, and denote by A/" the dimension of Z. We 
emphasize that the dimension Af is typically very large. These "truth" approximation 
subspaces inherit the inner products and norms of the exact spaces: (•, ■)x = (•, ■)x , 
II . \\x ^ II ■ ||;.„, (■, ■)y ^ i; ■)y., II • ||y ^ || • ||r., and (•, .)z ^ (•, ■)z., \\ ■ \\z ^ \\ ■ lU- 

Clearly, the continuity and coercivity properties (jl.ip . (|1.2p . and ()1.3p are passed 
on to the "truth" approximation spaces, 

7a(Ai) = sup sup < oo, y fieV, (1.6) 

uexvex ||m||x||w|U 

IbW = sup sup <oo, V ^ e X>, (1.7) 

qeYvex llglli-llwlU 

aaifi) ^ inf ^^^^ > 0, V /i e P, (1.8) 

vex \\v\\j^ 

and so is the inner product a(-, •; n) for any ^ E T>; thus, || • \\x,fi = \\ ■ \\x„,ti defines a 
norm on X that is equivalent to || • || . We now further assume that the approximation 
spaces X and Y are chosen such that they satisfy the Ladyzhenskaya-Babuska-Brczzi 
(LBB) inf-sup condition (see, e.g., [1]) 

(5Bri^i) = inf sup > PUf^) > 0, VfieV, (1.9) 

i^Yvex \\q\\Y\\v\\x 

where (/i) is a constant independent of the dimension Af. 

We now define our "truth" approximations to be the Galerkin projections of 
Uc(/i) G Xc and Po(/^) G onto X and Y, respectively: Given any /i G 2?, we find 
{u{fi) , p{fi)) e X xY such that 

a{u{n),v:fi) + b{v,p{fi);n)=f{v;fi), \/ v G X, 
b{u{n),q;n) = g{q;fi), V q e Y. 
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As for the exact problem in gTU it follows from ([L^ . ([TT7)) . ([TTS)) . and ([Q]) that the 
"truth" problem (|1.10p has a unique solution for any f{-;fi) E X^, .g(-;^) G FJ. The 
bilinear forms a{-, ■; /u) and •; /x) define bounded linear operators A{fi) : X — )> X' , 
B{ii) : X ^Y' and its transpose : F -> X' by 

u, v) ~ a{u, v; /i), y u,v E X, 

{B{^i)v,q)^h{v,q■^l)^{B{^lyq,v), yvex,yqeY; 

here, (■, ■} denotes the respective dual pairing. The "truth" system ()f .lOp can thus be 
equivalently written as 

A{^i)u{^i) + B{^lYp{^^) = f{^i) inx', 

B{^l)u{^l) = g{p) iny', 

where /(/i) = f{-\y)\x e X' and g{y) = g{-,y)\Y e Y' for all ^ e P. 

2. The Reduced Basis Method. We now turn to the reduced basis (RB) 
method. We here focus on the development of a priori and a posteriori error es- 
timates, which shall be specialized to the symmetric context. Other parts of the 
methodology such as the Offline-Online computational strategy as well as the con- 
struction of effective reduced basis approximation spaces shall only be briefly recalled 
as they have already been extensively discussed in [7]. 

2.1. Galerkin Projection. Suppose that we are given a set of nested, low- 
dimensional reduced basis approximation subspaces X^ C Xtv+i C X and Y^ C 
Yn+i c r, TV e n max = {Ij ■ ■ • ; -^max}- We denote by Nx and Ny the dimensions 
of X^r and Y/v, respectively, and the total dimension of = Xiq x Y/v by Nz = 
Nx+Ny- The subspaces Xx, Y/v, and Zx again inherit all inner products and norms 
of X, Y, and Z, respectively. The reduced basis approximation is then defined as the 
Galerkin projection onto these low-dimensional subspaces: For any given /x G 2?, we 
find MAr(/i) e Xx and pjv(m) G Yx such that 



(2.1) 



a(uAr(/i),UAr;/i) + b{vx ,PN{tj)\ = li'^N;^), V vx e Xx, 

b{uxi^J.),qx;^J■) = giqN^,^J'), ^qx^Yx. 

Written in operator notation, the discrete reduced basis system reads 

Ax{^JL)ux{^i)+ Bx{^if pxip) ^ fNUi) inX^, (2.2) 
Bx{tj)ux{y) = gx{y) inFj^,, (2.3) 

where fx{lj)=f{-]li)\xN G X^, ^^(m) = m)^™ S Fj^, and the bounded linear 
operators Ax{^j) ■ Xx — ^ X'j^, Bx{^j) ■ Xx Y^ and its transpose Bxif^Y ■ Yx — >■ 
X'j^ are given by 

{Ax{^J.)ux,vx) = a{ux,vx■■,^J,), W ux,vx G Xx, 

{Bxi^J,)vx,qx) = b{vx,qx;i-i') = {Bxi^J.Y qx,vx), ^ vx e Xx, V G Yx- 

We recall (see [7]) that the system (|2.1I) is well-posed if and only if the reduced basis 
approximation spaces Xx , Yx satisfy the inf-sup condition 

/3^(m)^ inf sup ii^^^'T'? >Q' (2.4) 

In&Yn ^^^Xn mN\\Y\\Vx\\X 

in this case, the pair (Xx,Yx) is called stable. 
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2.2. Reduced Basis Error Estimation. We here consider a priori as well as 
a posteriori estimates for the errors in the reduced basis approximations. 

In this section, we assume that the low-dimensional reduced basis spaces X^r, Y^r 
are constructed such that for any given parameter /i S I?, a solution (ujv(^),pAr(/i)) G 
Xn X Ym to (|2.1I) exists (see [3 §2.3]). We then denote the errors in the reduced basis 
approximations U7v(/^) G Xn, PwCm) ^ ^^v, and {uN{fJ-),PNitJ-)) G Zn with respect to 
the truth approximations by 

e'^iiJ.) = u{n) - UNifJ.) e X, 

eP^{^)=p{^)-pN{^) eY, (2.5) 
ejv(M) = (e]V(M),e^(A*)) e Z. 

2.2.1. j4 Priori Error Estimates. Concerning the rate at which the reduced 
basis approximations um{p-) and pn{p^) converge to the truth approximations u(/i) 
and p(/i), respectively, we can derive the following result. 

Proposition 2.1. For any given ^ £ V and N e Nmax, we have 



\\eUmx<2j^ inf \\u{^)-y^\\^ + 2!lil^ ini \\p{p) - q^lW; 

Sjv(^i)'^JV=ffJv(A') 



moreover, if the spaces Xpf^Y^ are stable (see ^2.1]) . we also obtain 



(2.6) 




qN&YN 



inf \\ui^i)-VN\\x. (2.7) 

-Bjv(A')i']v=ffjv(M) 

Proof. We here use techniques very similar to those presented in [4] for finite 
element approximations. Take any parameter fi G V and N G Nmax. Recall that 
a(-, •; fj.) defines an inner product on X and the associated norm is denoted by || • ||x,/i- 

By the definition of the reduced basis approximation in H2.1\ as the Galerkin 
projection of {u{fi) , p{ii)) onto Xn x Yn, the errors e^{fi) and e^(/i) satisfy 

a{e'^{fi),VN;fJ-) + b{vN,e''j^{p);tJ-) =0, VvNeXN- (2.8) 

First, we prove that p. 61) holds true. For any vn G -'^a' such that _Bjv(/i)vAr = 
gN{p^) in Y^, we have vn — un{p) G ker(i?Ar(/i)) and 

II .Ml a{vN -UNi^J-),VN -U]s[in);iJ,) 

\\VN - UN{mx,f. = r. 

\\VN - UNWWx.fi 

a{vN -UN{^l),vuN■,^J.) 
- n n 

tujvGkor(B„(/i)) \\WN\\X,tJ. 
t«jvekor(B„(/i)) llli'Afllx.At 

_ a{vN -u{^i),WN]tj) -b{wN,e%{iJL)\ij) 



t«]v6kcr(Bjv(M)) ||wAf||x,p 
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where the last equahty fohows from (|2.8I) . For wn G ker(i?jv(/i)), 6(i(jjv,pAr(/x); /i) = 
b{w]^, qn; fi) =0 holds for all qj^ G Y^r. Inserting this in the inequality above yields, 
for any qn G >Ar, 

\\VN -UN{mX,f, < sup r r 

u)ivekcr(S«(M)) ||W'A'|k,M 

< - u{^,)\\x., + ^li^lbOi) - gATllr, (2.9) 

where the latter is obtained from the Cauchy-Schwarz inequality for the inner product 
a(-, •; (|1.7p . and (jl.Sp . Using the triangle inequality and p.9p . 

l|e]V(Ai)l|x,M < [[-"(A*) - WAr||x,M + \\vN - ujv(a*)||x,;. 

< 2\\u{^,) - v^Wx,, + - qN\\Y, (2.10) 

the a priori stability estimate (12.61) then follows from (|2.10p . (|1.6p . and (|1.8p . 

We now turn to (|2.7p . Assuming that X^r and Y/v are stable, the inf-sup condition 
(|2.4p provides 

o / Ml / Ml ^ b{vN,qN -pN{fi);fi) ^ 

PAf(M)ll'?iv -pw(M)lk < sup 7 , VgN^YN. (2.11) 

vn&Xn W^nWx 

For any vn & and ^at € Y;v, we moreover have 

h{vN,qN -pNitA^iA = b{vN,qN -p(/i);/i) +b{vN,e^j^{p); p.) 

= b{vN, qN - p{p); fJ.) - a{e'%{n),VN; f^), 

where the last equality follows from (12. 8p . Applying this to (j2.1ip . together with the 
Cauchy-Schwarz inequality for a(-, •; /i) and (II. 6p . we obtain 



||g^-pA,(M)||y<^^||p(/.)-gA.|!y + ^^||e]tr(/i)|U,,„ VgArGyAr; (2.12) 

the a priori error estimate (j2.7p thus holds again from the triangle inequality, (|2.12p . 
dmni), and (HH). □ 

2.2.2. A Posteriori Error Estimates. We now show that the rigorous and 
computationally efficient reduced basis a posteriori error bounds A^^'^(^), A^^'^(/i), 
and A^''(/i) presented in [7] may be further sharpened in the case of a symmetric 
problem. 

Ingredients. To formulate rigorous and inexpensive upper bounds for the respec- 
tive errors defined in (12. 5p . we again rely on two sets of ingredients. The first set of 
ingredients consists of computationally (Online-)efficient lower and upper bounds to 
the truth continuity, coercivity, and inf-sup constants (|1.6I) . (jl.Sp . and (|1.9I) . 



^a^(M) < aa(/i) < 



a^:^(M) < «a(/i) < ar(A^), VmGI?. (2.13) 
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The second set of ingredients consists of dual norms of the residuals associated 
with the reduced basis approximation, 

||r^(-;A^)|U' = sup^i^, ^)||^_ sup ^i^^, (2.14) 

vex \\v\\x qeY \\q\\Y 

where, for all fJ. EV, rl^{-; ^j) G X' and r^(-; /i) G Y' are defined as 

rl^{v;n) = f{v;fi) ~ a{uN{fJ.),v;fi) ~b{v,pN{fJ.);f^), V v e X, (2.15) 

rNiiu^) ^ 9{q;f^) -b{'UN{f^),q;tJ'), VqeY. (2.16) 

Error Bound Formulation. We may now formulate a posteriori error bounds for 
the respective errors (j2.5p in the reduced basis approximation. 

Proposition 2.2. For any given p eV, N e Nmax, anda^^{p), Ja^ifJ-), /^ef (m) 
satisfying i2.13\) . we define 




A-»(.,.^^ + y^^|ip^, (2.17, 



Then, A^'''^™(/i) and A'^^^™(/i) are upper bounds for the errors e]^(/i) and e^(/i) 

||e5V(/.)||x < A™(^) < A^^^(^), ||e^^(^)||^ < < ^'/'M, (2-19) 

for all fi E T) and N G Nmax • 

Proof. Let /i be an arbitrary but fixed parameter in V and N G Nmax- We here 
proceed as in the proof of [71 Corollary 2.2], only that we may now exploit the fact 
that a(-, •; /i) defines an inner product on X; again, recall that the associated norm is 
denoted by || ■ \\ x,^l- 

By (PTT51) . (pflB)) . and prTU| . the errors e]^(/i) G X and e^(/i) G Y satisfy the 
equations 

a(e]^(/.i),t;;/i)+&(«,e^(M);Ai)=^^(«;A'), V « G X, (2.20) 
5(e5V(/^),g;A*)-r-2,(g;^), VgGF. (2.21) 

The error e^{fj,) G X in the approximation of the primal variable may now be 
uniquely decomposed into e'%{y) = e%{fi) + ej^ifi) where e%{fi) G ker(i?(/i)) and 
ej^il^) G X such that 

a{e^{fi),vo;fi) = 0, Vt^o Gker(S(Ai)). (2.22) 
From (p:^ and e%{^i) G ker(B(^)) then solves 

a{e%{fi),vo; = r]v(t;o;/i) - a(e]^(/i), uq; m) = ''7v(^^o;m), V G ker(B(/i)). 
Setting here = e^(/i), we have 

P°A^(M)ll^,^-rj,(g°^(M);/i) < We'MWx sup ^i^^ 
< / — ^ l|ejv(wlU,M sup — — , 

Vaa(M) t>oGkor(S(M)) IPo||X 
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where the last inequahty follows from (ll.Sp . Hence, e'j^in) is bounded by 



|e^(MjIU,,. < — ^== sup —rr—r < 



To obtain an upper bound for ej^{fi), we here consider the inf-sup constant 

;9(M)^inf sup . (2.24) 

leYyfzx \\q\\Y\\v\\x,t^ 



From (11.6P and (ll.Sp . we have 



^^^^''^ </3(m)<^^, (2.25) 



and thus particularly /3(/i) > by the LBB inf-sup condition (|1.9p . Consequently, 

P{mv\\x,t. < sup— — 
q&Y \\q\\Y 

holds true for any v G X such that a{v, vq; h) = for all vq £ ker(_B(^)) (see, e.g., [U 
§11.1, Proposition 1.2]). Applied to ej::i{ij) satisfying (|2.22p . this yields 

ll-;-L/,^ll ^ 1 „„„ Kgj^(M),g;A^) _ lk^(-;M)l|y' 

e^lA*) < sup ^— r ' (2.26) 

where the equality follows from (12.211) as B{^) ej^{n) = -B(/i) e'^{fi) in y . Now, using 
the triangle inequality, we may finally derive that 

~o /.,\ 



\\e%{^^)\\x,,<\\iM\U, + \\eMU 



by combining ((2^, (|2:26)) . and ([223; the bound (|2T7)) thus follows from (EJH), 
(OD . and (EISl). 

For the error e^(/i) in the approximation of the Lagrange multiplier, we obtain 
from dm]) and IJ^L^ that 



-N\l^)\\y — '3 TT sup TT—r. — — 1-^ sup - - 

PbtW vex \\v\\x PbAn vex \\v\\x 



- ^^("'^^^■'''^"''' ^ (2.28) 

where the last inequality holds by the Cauchy-Schwarz inequality for the inner product 
a(-, •; A^) and pT^ . Together with ^(TFJ) and again (I^T^ . this leads to (HjH]). 

As it is clearly aa(/x) < 7a (a*) from the definitions in (|1.6p and (|1.8p . we have 
A^'^y'"(Ai) < A];;^' (Ai) and A^^y'"(/x) < A^^'^(Ai); this eventually yields (EH]). □ 

Remark 2.3. We here comment on the above error bounds in a little more detail. 

First, we clearly obtain 



A?^-(m) ^ v^(A™(a.))' + (A^^^-(^))^ ^eP, iVeN^a., (2.29) 
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as a rigorous upper bound for the combined error €^{(1) such that 

\\eN{pi)\\z < A5^'"(/i) < AB'(Ai), V A* e 2?, V TV e N„,ax. 
Analogous to [71 Remark 2.3], the associated efFectivities satisfy 

1 < ^JT\^,) = y < C(A*)(7a(A*) + Ibifi)), V e P, V TV e N„,ax, (2.30) 

\\eNW\\z 

where C(^) > is a constant depending on aj;^(/z), 7y^(At), and /^ef (m)- 

Second, also the error bound (2.22) given in [71 Remark 2.3] may be specialized 
to the symmetric case. Let /i be any parameter in V. Using p.23p . ()2.26p . ()2.25p . 
(fLSl) . and 

||r^(-;M)Ilx'<7a(A^)l|e]^(M)llx+76(M)l|e^(M)||„ 
\K{■■,^^)\\Y' <lbif^)\\eUmx, 

we obtain that 



f]v(^o;Ai) , /7a(Ai) lkw(-;M)lh 



aa(Ai) i;oekcr(B(Ai)) ll«o|U V "a (m) /^Br (m) 

represents a reliable and efficient a posteriori error estimator only for ||e5v(/^)ll^- 

Third, the error e]^(/^) may clearly also be measured in the energy norm || • \\x,n — 
Va{-, SM)- For any given fi £V, N e N^ax, and a^^{n), -f™{fJ.), ^ef (m) satisfying 
(Pl^ . we define 



X«,sym/ \_ lkA'('j M) I yla^ ifJ-) n 2/ Ml /o 01 \ 

A. (^) = ^7^m^ + 7LB(^IM-;^)ll- (2.31) 



Then, from (j2.27p and (|2.13p . A'^''^™{fi) is an upper bound for e^lji) such that 

||e]V(M)llA> < A;;;^y'°(/i), e 2?, VTV e N„,ax. (2.32) 



2.3. Offline— Online Computational Procedure. The efficiency of the re- 
duced basis method relies on an Offline-Online computational decomposition strat- 
egy. As this is by now standard, we shall only provide a brief summary at this point 
and refer to, e.g., |14) for further details. The procedure requires that all involved op- 
erators can be afRnely expanded with respect to the parameter fi. All /i-independent 
quantities are then formed and stored within a computationally expensive Offline 
stage, which is performed only once and the cost of which depends on the large truth 
dimension TV. For any given fj, G V, the reduced basis approximation is then com- 
puted within a highly efhcient Online stage; the cost does not depend on TV but only 
on the considerably smaller dimension of the reduced basis approximation space. 

The computation of the a posteriori error bounds clearly consists of two com- 
ponents: The calculation of the residual dual norms (I2.14p and the lower and upper 
bounds (|2.13p to the involved coercivity, continuity, and inf-sup stability constants. 
The former is again an application of standard reduced basis techniques that can 
be found in [T3]. The latter is achieved by a successive constraint method (SCM) 
proposed in [S]; we refer to [6j for details in our saddle point context. 
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2.4. Construction of Reduced Basis Approximation Spaces. The low- 
dimensional reduced basis approximation spaces X^Vj^at, N g Nmax, are constructed 
by exploiting the parametric structure of the problem: According to the so-called 
Lagrange approach, basis functions are given by truth solutions associated with sev- 
eral chosen parameter snapshots. In the case of saddle point problems, additional 
care must be taken in the construction of stable reduced basis approximation spaces 
(see §2.1|) . This has been extensively discussed in [7]: Stability is achieved through 
enriching the reduced basis space for the primal variable appropriately. Different 
strategies may be applied: We may add supremizer functions [131 115] or additional 
truth solutions ^ favoring either the approximations for the primal or the Lagrange 
multiplier variables, respectively. 

Keeping computational cost to a minimum, we aim to construct stable approxi- 
mation spaces Xn^Yn that appropriately represent the submanifold associated with 
the parametric dependence with as few basis functions as possible. To this end, we 
invoke a greedy sampling process [H [S] in which our rigorous and computationally 
inexpensive reduced basis a posteriori error bounds are used to identify truth solu- 
tions that are not yet well-approximated. In [7], this procedure has been extended to 
the special needs of saddle point problems. Correspondingly, we here again consider 
the following three sampling processes: the standard greedy sampling process sum- 
marized in [71 Algorithm 1] as well as the modified procedures [71 Algorithm 2] and 
[3 Algorithm 3] where the need for stabilization is recognized adaptively. 

3. Numerical Results. We now apply the reduced basis methodology devel- 
oped in ij2] to the Stokes model problem described in [6l [7] . Numerical results are 
attained using the open source software rbOOmit [ir , an implementation of the re- 
duced basis framework within the C-f-|- parallel finite element library libMesh \T0*. 

Motivated by applications in the field of microfluidics, we consider a Stokes flow in 
a two-dimensional microchannel with a parametrized, rectangular obstacle (see [7]). 
As our truth discretization in §1.21 '^6 then choose a mixed finite element method 
using the standard conforming P2-P1 Taylor-Hood approximation spaces [16^; the 
truth system (|1.10p here exhibits a dimension oi J\f = 72,076. 

Since our Stokes model problem is clearly symmetric, we may now compare the re- 
duced basis a posteriori error bounds A'^^^'^{^), A^^^^"^{fi), and A^™(^) (see (|2.17p . 
()2.18|) . and (|2.29p ) specific to the symmetric case with the more general bounds 
A]y "^(/i), "^(m); and A^''(/i) introduced in 7, Corollary 2.2]; in addition, we shall 
again consider the reduced basis error estimate A^^{fi) based on Babuska's theory for 
noncoercive problems as defined in [71 Corollary 2.1]. For this purpose, we build the 
approximation spaces Xj^ and Ypf by using either Algorithm 1, Algorithm 2, or Algo- 
rithm 3 (see > j2.4[) . All sampling procedures are based on an exhaustive random sample 
S cPofsizelSj =4,900 and the relative error bound Ajv(a*) = ^^''^™{fJ-)/\\uN{fJ-)\\x 
(see (|2.31|) ): in Algorithm 2 and Algorithm 3, we further set S^^^ — 0.1. 

Figure 13.11 and Figure 13.21 show the maximum errors in the reduced basis ap- 
proximations for the primal variable and the Lagrange multiplier, respectively, to- 
gether with the respective error bounds A^^^'^di), A^^'^(/i), A^''^™(^), A^j^^'^{fi), 
and A^^(^); Fi gure 13.31 shows the maximum total error in the reduced basis approx- 
imation and associated error bounds A^™(/i), A^'^(/i), and A^^(/i). The maximum 
is computed over a sample of 25 parameter values. For this sample, to analyze only 
the effects of the a posteriori error bound formulations and eliminate contributions 
of the SCM (see > j2.3|) . we again use the exact constants rather than the lower/upper 
bounds (|2.13p . Effectivities associated with the error bounds are given in Table 13.11 
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(a) (b) (c) 

Fig. 3.1. Maximum relative error \\e'^{fj,)\\x /\\'a{iJ,)\\x '""^ maximum relative error bounds 
A^''^™(At)/||n{/i)||x, A'^^''{fj.)/\\u{ii)\\x , and A^''{ij,)/\\u{^i)\\x shown as functions of Nz for 
(a) Algorithm 1, (b) Algorithm 2, and (c) Algorithm 3 (see %2.4'i ; the maximum is taken over 
25 parameter values; the computation of the error bounds is based on the exact constants involved. 




(a) (b) (c) 

Fig. 3.2. Maximum relative error \\e^{fi)\\Y /\\pitJ')\\Y cii^d maximum relative error bounds 
A^''^™{/i)/||p{/i)||y, A^^''(/i)/||p(^)||y , and Aff'{fi)/\\p(ji)\\Y shown as functions of Nz for (a) Al- 
gorithm 1, (b) Algorithm 2, and (c) Algorithm 3 (see %2.Ji.^ ; the maximum is taken over 25 parameter 
values; the computation of the error bounds is based on the exact constants involved. 



As before (see [7]), maximum effectivities associated with "^(/i) essentially range 
from 50 to 80 (see Fig. 13.11 and Table 13. ip . Now, exploiting the symmetry of the 
problem, A^^'*'™(/x) provides a sharper bound not only in theory (see Proposition l2.2p 
but also in practice: Here, effectivities essentially vary between 20 and 40. Further- 
more, it is not surprising to observe that the best results are achieved by the bound 
^ti^bym^^^ that overestimates the error in the reduced basis approximations for the 
primal variable measured in the energy norm (see (12.321) ) only by a factor of approx- 
imately 15 (see Table lO]) . For A^'''*'™(/i), we still obtain rather large effectivities 
with average values of 362, 161, and 127 in case of Algorithm 1, Algorithm 2, and 
Algorithm 3, respectively; however, compared to A^^'(/x), the improvement is sig- 
nificant (see Fig. 13.21 and Table [5T|) . We emphasize that A^''-*""(^) and A^'^''"^(^) 
do in neither case perform worse than A^^(/i) but generally provide bounds for the 
errors in the primal and Lagrange multiplier variables that are much more accurate 
(see Table EH). 

We now discuss the Online computation times of the proposed methods. The 
SCM (see i i2.3|) enables the (Online-) efficient estimation of the constants aa(/i), 7a (m), 
and /^Br (m) ■ We here apply the method with the configurations as specified in and 
receive very accurate (Online-)efficient bounds Q^{p), 7y^(^), and providing 
a posteriori error bounds A]^(/x), A^''^™(/i), A^''^™(/i) and A^(/x), A5Jj-*''^™(/x) that 
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(a) (b) (c) 

Fig. 3.3. Maximum relative error \\e]N{fi)\\z /\\{u{fJ.),p{fi))\\z cin-d maximum relative error 
bounds A^'"(/i)/||(n{/i),p(At))||z, Afj" {^l)/\\{u{^l), plij,))\\z , and A^''{fj.)/\\{u{n), p{^i))\\z shown as 
functions of Nz for (a) Algorithm 1, (b) Algorithm 2, and (c) Algorithm 3 (see j jg.^l ); the maxi- 
mum is taken over 25 parameter values; the computation of the error bounds is based on the exact 
constants involved. 



(a) Algorithm 1 



N 


Nz 


~if,sym 

Vn 


ii,sym 

Vn 




Vn 


„P,sym 

Vn 


V% 


^« 


sym 

v^ 


„Br 

Vn 


Ba 

Vn 


5 


15 


13.21 


20.68 


58.04 


353.6 


202.1 


1534 


296.4 


148.0 


1238 


203.6 


10 


30 


12.34 


19.94 


63.93 


277.5 


292.5 


2268 


326.5 


168.6 


1559 


186.7 


15 


45 


13.22 


19.57 


57.42 


309.9 


331.5 


3223 


417.2 


201.2 


1598 


228.6 


20 


60 


13.29 


19.35 


58.40 


255.0 


355.9 


3240 


415.7 


208.1 


1580 


207.1 


25 


75 


11.97 


18.77 


57.92 


242.9 


405.8 


3513 


476.3 


198.4 


1560 


216.4 


30 


90 


13.55 


19.93 


54.23 


282.3 


399.2 


3423 


453.6 


227.0 


1678 


232.6 


35 


105 


13.62 


20.44 


53.93 


281.7 


432.1 


2502 


519.8 


235.2 


1535 


234.9 


40 


120 


14.28 


21.53 


57.49 


295.4 


482.6 


3317 


495.2 


253.0 


1526 


253.1 



(b) Algorithm 2 



N 


Nz 


-u.sym 

Vn 


Ti .sym 

Vn 


U.Br 

Vn 


u.JrJa 

Vn 


„P.sym 

Vn 


V% 


p.JrJa 

Vn 


sym 

Vn 


„Br 

Vn 


Ba 

Vn 


5 


12 


16.11 


23.89 


55.18 


715.1 


98.47 


756.3 


142.5 


88.04 


673.4 


139.7 


10 


23 


22.45 


32.53 


62.65 


1894 


121.8 


799.9 


234.8 


99.72 


700.7 


233.0 


15 


34 


13.57 


25.31 


68.15 


425.1 


152.8 


766.6 


253.9 


130.7 


660.6 


216.6 


20 


46 


13.17 


20.31 


61.41 


505.3 


224.0 


1268 


307.0 


172.7 


1081 


236.1 


25 


56 


16.44 


27.51 


63.04 


1013 


97.71 


550.1 


258.9 


94.87 


523.1 


250.8 


30 


68 


14.63 


23.48 


61.14 


583.5 


237.9 


1287 


294.2 


189.6 


1023 


233.8 


35 


78 


14.80 


22.35 


59.63 


707.3 


126.3 


843.6 


269.0 


118.3 


778.2 


251.5 


40 


89 


14.13 


23.70 


61.52 


613.1 


162.4 


800.6 


285.4 


139.3 


685.6 


244.3 


(c) Algorithm 3 


N 


Nz 


-u.sym 

Vn 


Ti .sym 

Vn 


U.Br 

Vn 


li.JrJa 

Vn 


„P.sym 

Vn 


p,Br 
'"fN 


p.JrSa 

v'n 


sym 

vZ 


„Br 

Vn 




5 


13 


18.36 


32.56 


68.23 


1146 


85.38 


793.3 


162.2 


79.09 


755.4 


157.9 


10 


29 


22.08 


33.22 


72.73 


717.1 


183.0 


996.8 


219.0 


152.9 


830.9 


182.5 


15 


39 


15.56 


26.58 


63.64 


446.0 


137.4 


721.9 


195.0 


123.8 


649.1 


173.1 


20 


49 


15.13 


24.88 


69.96 


699.3 


118.1 


807.8 


174.9 


104.2 


753.7 


162.0 


25 


60 


19.45 


33.86 


71.13 


1034 


130.3 


577.4 


206.5 


115.1 


509.1 


195.4 


30 


73 


20.00 


40.18 


97.51 


976.6 


102.5 


582.7 


187.7 


94.86 


538.1 


179.6 


35 


83 


16.90 


28.72 


71.93 


821.5 


136.3 


692.0 


217.4 


124.6 


631.3 


198.3 


40 


93 


17.87 


33.31 


82.45 


973.4 


196.7 


1012 


301.9 


165.4 


849.5 


288.4 



Table 3.1 

Maximum effectivities rj^'''^™(/i) = A^''^™(/i)/|le]^(/i)||x,^ (see 112.32( 1). Vn'^^'^M = 
A^^^y'^i^^)/\\eli^,)\\x, V^N^'^'^M = A^^^y^i^,)/\\eW\\Y, v'TM = A=^-(/.)/||e^(M)||z (see 
mB and ^EM), an<ir,]^'S'(/.) = A"^''\^^)/\\e^^i^,)\\x , rj^^^^M = A'}}'"^ i^^)/\\e%i^,)\\Y , v'^'M = 
A^^i^^)/\\eN(p)\\z, V^^^'^if^) = A%-(p)/\\el(p)\\x, = A^-M/\\eW\\Y , vW(P) = 

^^^M/\\<^N{tJ-)\\z (see ^) for (a) Algorithm 1, (b) Algorithm 2, and (c) Algorithm 3 (see ^2.4^ ; 
the maximum is taken over 25 parameter values; the computation of the error bounds is based on 
the exact constants involved. 
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essentially coincide with their values based on the evaluation of the exact constants; 
associated effectivities thus remain the same as shown in Table 13.11 For comparison, 
once the /i-independent parts in the affine expansion of the involved operators (see 
i|2.3p have been formed, direct computation of the truth approximation {u{fi),p{n)) 
(i.e., assembly and solution of (ILIOI) ) takes on average 6.5 seconds on a 2.66 GHz Intel 
Core 2 Duo processor. The rigorous and (Online-) efficient error bounds A^^^"^(fi) and 
^Pj^symji^'j allow us to choosc the reduced basis system dimension Nz just large enough 
to obtain a desired accuracy. In case of Algorithm 1, we need Nz = 51 to achieve 
a prescribed accuracy of roughly 1% or better in the reduced basis approximations 
UAr(/i) (see Fig. 13. U fa)). Once the database has been loaded, the Online calculation of 
{uff{^),Pff{^)) (i.e., assembly and solution of (12. ip ) and A'^^^'^{fi), A^''^™(/i) for any 
new value of /i G P takes on average 0.31 and 20.99 milliseconds, respectively, which is 
in total roughly 300 x faster than direct computation of the truth approximation. We 
note that again (see 7 ), stabilizing adaptively pays off: In the case of Algorithm 3, 
the same accuracy is achieved for Nz = 35 (see Fig. 13.11 (c)); the Online calculation 
of (wAr(/x),pAr(/i)) and A^^^'^{fi), A^jf^'^lfi) then takes on average 0.14 and 13.29 
milliseconds, respectively, and is thus roughly 480 x faster than direct computation 
of the truth approximation. Detailed computation times, also for Algorithm 2, are 
given in Table 13.21 Clearly, due to the significant improvement of the new error 
bounds in terms of sharpness, the new methods guarantee a prescribed accuracy in 
the reduced basis approximations at notable Online savings when compared to the 
methods presented in [7] (see Table and [71 Table 2.3]). 



(a) Accuracy of at least 1% (resp., 0.1%) for the reduced basis approximations Uff(fj,) 



Method 


Nz 


N 


("Jv(M).P]v(/i)) 


'^N (m) 


Total 


Algorithm 1 


51 (78) 


17 (26) 


0.31 (0.79) 


20.99 (39.60) 


21.20 (40.38) 


Algorithm 2 


44 (61) 


19 (27) 


0.21 (0.42) 


16.25 (24.73) 


16.46 (25.15) 


Algorithm 3 


35 (55) 


13 (23) 


0.14 (0.32) 


13.29 (21.59) 


13.43 (21.92) 


(b) Accuracy of at least 1% (resp., 0.1%) for the reduced basis approximations PjvCm) 


Method 


Nz 


TV 


("Jv(M).P]v(/i)) 




Total 


Algorithm 1 


54 (84) 


18 (28) 


0.35 (0.91) 


22.75 (44.70) 


23.10 (45.61) 


Algorithm 2 


44 (72) 


19 (32) 


0.21 (0.58) 


16.25 (31.54) 


16.46 (32.12) 


Algorithm 3 


37 (62) 


14 (26) 


0.16 (0.44) 


13.94 (25.77) 


14.09 (26.21) 



Table 3.2 

Average computation times in milliseconds for the Online evaluation of {ui^{fi),pisr{fj,)) (assem- 
bly and solution of i2.1\ }) and the error bounds A^'^^'^{fi), A^'''*'™(/i) (see \2.17^ , i2. 18V ): times 
are measured using either Algorithm 1, Algorithm 2, or Algorithm 3 (see with a prescribed 

accuracy of 1% (resp., 0.1%) for the reduced basis approximations (a) uj^{^) and (b) ppf^fi). 



4. Conclusion. In this paper, we improve reduced basis error bounds introduced 
in [7] for the special case of a symmetric problem. Numerical results provide a direct 
comparison with former approaches and demonstrate the superiority of the developed 
bounds with respect to sharpness; in particular, the upper bounds provided for the 
errors in the reduced basis approximations for the primal variable exhibit effectivities 
that are comparatively low. As a direct consequence, the proposed methods provide 
accurate reduced basis approximations at much less computational cost. 
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